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Abstract 



Both the mass spectra and the wave functions of the light pseudoscalar 
(n, K, 77, 7/) and vector(p, K*,u, </>) mesons are analyzed within the framework 
of the light-cone constituent quark model. A gaussian radial wave function 
is used as a trial function of the variational principle for a QCD motivated 
Hamiltonian which includes not only the Coulomb plus confining potential 
but also the hyperfine interaction to obtain the correct p — tt splitting. For 
the confining potential, we use (1) harmonic oscillator potential and (2) linear 
potential and compare the numerical results for these two cases. The mixing 
angles of cu — (ft and r] — rj' are predicted and various physical observables such 
as decay constants, charge radii, and radiative decay rates etc. are calculated. 
Our numerical results in two cases (1) and (2) are overall not much different 
from each other and have a good agreement with the available experimental 
data. 
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I. INTRODUCTION 



It has been realized that the relativistic effects are crucial to describe the low-lying 
hadrons made of u,d and s quarks and anti-quarks The light-cone quark model [@-[12| 
takes the advantages of the equal light-cone time(r = t + z/c) quantization and includes 
the important relativistic effects in the hadronic wave functions. The distinct features of 
the light-cone equal-r quantization compared to the ordinary equal-t quantization may be 
summarized as the suppression of vacuum fluctuations with the decoupling of complicated 
zero-modes and the conversion of the dynamical problem from boost to rotation. 

The suppression of vacuum fluctuations is due to the rational energy-momentum dis- 
persion relation which correlates the signs of the light-cone energy k~ = k° — k 3 and the 
light-cone momentum k + = k° + k 3 ||. However, the non-trivial vacuum phenomena can still 
be realized in the light-cone quantization approach if one takes into account the non-trivial 
zero-mode(A; + = 0) contributions. As an example, it is shown that the axial anomaly in 
the Schwinger model can be obtained in the light-cone quantization approach by carefully 
analyzing the contributions from zero-modes. Therefore, in the light-cone quantization ap- 
proach, one can take advantage of the rational energy-momentum dispersion relation and 
build a clean Fock state expansion of hadronic wave functions based on a simple vacuum by 
decoupling the complicated non-trivial zero-modes. The decoupling of zero-modes can be 
achieved in the light-cone quark model since the constituent quark and anti-quark acquire 
appreciable constituent masses. Furthermore, the recent lattice QCD results |T3| indicated 
that the mass difference between rj and pseudoscalar octet mesons due to the complicated 
nontrivial vacuum effect increases (or decreases) as the quark mass m q decreases (or increases), 
i.e., the effect of the topological charge contribution should be small as m q increases. This 
supports us to build the constituent quark model in the light-cone quantization approach 
because the complicated nontrivial vacuum effect in QCD can be traded off by the rather 
large constituent quark masses. One can also provide a well-established formulation of var- 
ious form factor calculations in the light-cone quantization method using the well-known 
Drell-Yan-West(g + = 0) frame. We take this as a distinctive advantage of the light-cone 
quark model. 

The conversion of the dynamical problem from boost to rotation can also be regarded 
as an advantage because the rotation is compact, i.e., closed and periodic. The reason why 
the rotation is a dynamical problem in the light-cone quantization approach is because the 
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quantization surface r = is not invariant under the transverse rotation whose direction is 



perpendicular to the direction of the quantization axis z at equal r [To] . Thus, the transverse 
angular momentum operator involves the interaction that changes the particle number and 
it is not easy to specify the total angular momentum of a particular hadronic state. Also r 
is not invariant under parity |IE . We circumvent these problems of assigning the quantum 
numbers J PC to hadrons by using the Melosh transformation of each constituents from equal 
t to equal r. 

In our light-cone quark model of mesons, the meson state \M > is thus represented by 

\M>=^ Q \QQ>, (1.1) 

where Q and Q are the effective dressed quark and anti-quark. The model wave function is 
given by 



dk 

= *{x, k ± , A„ Xq) = y -£4>{x, k ± )TZ(x, k ± , X q , A,), (1.2) 

where 0(x, kj_) is the radial wave function, dk n /dx is a Jacobi factor and 7Z(x, kj_, X q , Xq) is 
the spin-orbit wave function obtained by the interaction-independent Melosh transformation. 
When the longitudinal component k n is defined by k n = (x — l/2)Mo + (m| — m q )/2Mo, the 
Jacobian of the variable transformation {x, k±} — > k = (k n , kj_) is given by 



8k n M < v(m^-m^2 



dx 4x(l — x) 



1 



\™>q ~ m qJ 



M, 



2 



;i.3) 



The explicit spin-orbit wave function of definite spin (S, S z ) can be obtained by 
TZ(x,k ± , XgXq) = ^2 < Xq\lZl A (x,\i ± ,m q )\s q > 

x < Aq|7^i I (l - x, -k ± ,rriq)\sq >< ]-S q -Sq\SS z >, (1.4) 

where the Melosh transformation is given by 

m + xM -ia -(hxk) 
n M (x,k ± ,m) = (1.5) 

(m + xM Q ) 2 + k\ 



with n = (0, 0, 1) being a unit vector in the ^-direction. 

While the spin-orbit wave function is in principle uniquely determined by the Melosh 
transformation given by Eq.(1.5), a couple of different schemes for handling the meson mass 
Mq in Eq.(1.5) have appeared in the literatures 0HT^[. While in the invariant meson mass 
scheme |]2],|6Hl2l, the meson mass square Mq is given by 
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2 _ ki + m\ k 2 ± + m 



M = + (1.6) 

x 1 — x 



in the spin-averaged meson mass scheme [§-0], M was taken as the average of physical 
masses with appropriate weighting factors from the spin degrees of freedom. Nevertheless, 
once the best fit parameters were used |5]||, both schemes provided the predictions that 
were not only pretty similar with each other but also remarkably good |||7j compared to the 



available experimental data |17] for form factors, decay constants, charge radii etc. of various 
light pseudoscalar(7r, K, rj, rj') and vector(p, K*,ou, <fi) mesons as well as their radiative decay 
widths. The main difference in the best fit parameters was the constituent quark masses, 
i.e., m u = rrid= 330 MeV, m s = 450 MeV in the spin-averaged meson mass scheme |0-|5| 
while m u = m^= 250 MeV, m s = 370 MeV in the invariant meson mass scheme J7|. 

Also, among the literatures |7|,§,|nj using the invariant meson mass scheme, some litera- 
tures |0,H used the Jacobi factor dk n /dx in Eq.(1.2) while some did not. However, we 
have recently observed || that the numerical results of various physical observables from 
Refs. were almost equivalent to those of Ref. 0] regardless of the presence-absence of 
the Jacobi factor if the same form of radial wave function (e.g. Gaussian) was chosen and 
the best fit model parameters in the radial wave function were used. 

However, the effect from the difference in the choice of radial wave function, e.g., har- 



monic oscillator wave function plJS.lO] versus power-law wave function [11], was so substan- 



tial that one could not get the similar result by simply changing the model parameters in the 
chosen radial wave function. For example, in the phenomenology of various meson radiative 
decays at low Q 2 , we observed || that the gaussian type wave function was clearly better 
than the power-law wave function in comparison with the available experimental data. On 
the other hand, the radial function so far has been mostly taken as a model wave function 
rather than as a solution of QCD motivated dynamical equation. Even though the authors 
in Ref. [12] adopted the quark potential model developed by Godfrey and Isgur |IJ to repro- 
duce the meson mass spectra, their model predictions included neither the mixing angles of 
uj — cf) and 77 — 77' nor the form factors for various radiative decay processes of pseudoscalar 
and vector mesons. 

In this work, we are not taking exactly the same quark potential developed by Godfrey 
and Isgur [l| . However, we attempt to fill this gap between the model wave function and the 
QCD motivated potential, which includes not only the Coulomb plus confining potential but 
also the hyperfine interaction to obtain the correct p— ir splitting. For the confining potential, 
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we take (1) harmonic oscillator(HO) potential and (2) linear potential and compare the 
numerical results for these two cases. We use the variational principle to solve the equation 
of motion. Accordingly, our analysis covers the mass spectra of light pseudoscalar(7r, K, rj, rf) 
and vector(p, K*, u, <fi) mesons and the mixing angles of u — <p an d i] — rf as well as other 
observables such as charge radii, decay constants, radiative decay widths etc.. We exploit 
the invariant meson mass scheme in this model. We also adopt the parametrization to 



incorporate the quark-annihilation diagrams |18H20f mediated by gluon exchanges and the 
SU(3) symmetry breaking, i.e., m u ^) 7^ m s , in the determination of meson mixing angles. 

The paper is organized as follows: In Sec. II, we set up a simple QCD motivated effective 
Hamiltonian and use the gaussian radial wave function as a trial function of the variational 
principle. We find the optimum values of the model parameters, quark masses(m n (^, m s ) 
and gaussian parameters( (3 uU = [3 U ^ = P dd , j3 uS , /3 S g) for the two cases of confining potentials 
(1) and (2). We also analyze the meson mass spectra and predict the mixing angles of 
u — <j) and rj — rj'. We adopt the formulation to incorporate the quark- annihilation diagrams 
and the effect of SU(3) symmetry breaking in the meson mixing angles. In Sec. Ill, we 
calculate the decay constants, charge radii, form factors and radiative decay rates of various 
light pseudoscalar and vector mesons and discuss the numerical results of the two confining 
potentials (1) and (2) in comparison with the available experimental data. Summary and 
discussions follow in Sec. IV. The details of fixing the model parameters and the mixing angle 
formulations are presented in Appendices A and B, respectively. 

II. MODEL DESCRIPTION 

The QCD motivated effective Hamiltonian for the description of the meson mass spectra 
is given by 0,0 



n qq\^nlm > 



m 2 + k 2 + JrrA + k 2 + V ai 



I nlm i 



Hn + V 



>= M qq -\V s n t >, (2.1) 



where M qq - is the mass of the meson, k 2 = ~k\ + k\, and \^nim > * s the meson wave func- 
tion given in Eq.(1.2). In this work, we use the two interaction potentials V qq - for the 
pseudoscalar(0 h ) and vector(l ) mesons:(l) Coulomb plus harmonic oscillator (HO), and 
(2) Coulomb plus linear confining potentials. In addition, the hyperfine interaction, which 
is essential to distinguish vector from pseudoscalar mesons, is included for both cases, viz. 
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V qq = V (r) + V hyp (r) 

= a + V con , - ^ + ^W Coull (2.2) 

where V CO nf. = br[r 2 ] for the linear[HO] potential and < S q ■ Sq >= l/4[— 3/4] for vec- 
tor [pseudoscalar] meson. Even though more realistic solution of Eq.(2.1) can be obtained 
by expanding the radial function <fin,i=o{k 2 ) onto a truncated set of HO basis states |1,12 



? - e -j Y,n=T c n(fin,o(k 2 ) , our intention in this work is to explore only the 0~ + and 1 ground 
state meson properties. Therefore, we use the IS" state harmonic wave function lo (fc 2 ) as a 
trial function of the variational principle 

M*,k.i) = (-^) 1/2 exp(-/c72/3 2 ), (2.3) 

where <fi(x, kj_) is normalized according to 

X;^ 1 ^ y rf 2 k^|^f ^(a:, k^,zvz7)| 2 = ^'rfa; y rf 2 k_L | ^ 10 (a: , ) | 2 = 1. (2.4) 

Because of this rather simple trial function, our results could be regarded as crude approx- 
imations. However, we note that this choice is consistent with the light-cone quark model 
wave function which has been quite successful in describing various meson properties P~|TU 



Furthermore, Eq.(2.3) takes the same form as the ground state solution of the HO poten- 
tial even though it is not the exact solution for the linear potential case. As we show in 
Appendix A, after fixing the parameters a, b and k, the Coulomb plus HO potential Vo(r) 
in Eq.(2.2), turns out to be very similar in the relevant range of potential(r < 2 fm) to the 
Coulomb plus linear confining potentials [see Figs. 1(a) and 1(b)] which are frequently used 
in the literatures [|T1,|T2,21-25||. The details of fixing the parameters of our model, i.e., quark 
masses (m^), m s ), gaussian parameters(/3 M( j, j3 uS , j3 sS ) and potential parameters (a,b,n) in 
V q q given by Eq.(2.2) are summarized in the Appendix A. 

Following the procedure listed in the Appendix A, our optimized model parameters are 
given in Table I. In fixing all of these parameters, the variational principle [Eq. (A 1)] plays 
the crucial role for ud, us, and ss meson systems to share the same potential parameters 
(a,b,n) regardless of their quark- ant iquark contents[see Figs. 2(a) and 2(b)]. 

We also determine the mixing angles from the mass spectra of (u, <j)) and (r),rj'). Iden- 
tifying (fi, f 2 ) = ((fi,u)) and (r],r]') for vector and pseudoscalar nonets, the physical meson 
states fi and fi are given by 
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\fi > = — sin 5\nn > — cos S\ss >, 

|/ 2 > = cos5|nn > — sin 6\ss >, (2.5) 

where \nn >= l/y/2\uu + dd > and § = 9su(3) ~ 35.26° is the mixing angle. Taking into 
account SU(3) symmetry breaking and using the parametrization for the (mass) 2 matrix 
suggested by Scadron PD[ , we obtain 

tan2j (M%-M^(M%-Ml) 

ta 6 ~ {M%-Ml){Ml-Mt^ (2 ' 6) 

which is the model independent equation for any meson qq nonets. The details of obtaining 
meson mixing angles using quark-annihilation diagrams are summarized in the Appendix 
B. In order to predict the 10 — (j) and rj — rj' mixing angles, we use the physical masses [17] 



of Mfa = (m^m,) and Mj 2 = (m^rry) as well as the masses of M^ s = 996[952] MeV and 
Mf s = 732 [734] MeV obtained from the expectation value of H sS in Eq.(2.1) for HO [linear] 
potential case [see Appendix A for more details]. Our predictions for u — and 77 — 77' 
mixing angles for HOflinear] potential are \5y\ ~ 4.2°[7.8°] and 8su{3) ~ — 19.3°[— 19.6.°], 
respectively. The used mass spectra of light pseudoscalar and vector mesons are summarized 
in Table II. Since the signs of 5y for to — (p mixing are not yet definite 1181-^01,^71-1291 in the 
analysis of the quark-annihilation diagramfsee Appendix B], we will keep both signs of 5y 
when we compare various physical observables in the next section. 



III. APPLICATION 



In this section, we now use the optimum model parameters presented in the previous 
section and calculate various physical observables; (1) decay constants of light pseudoscalar 
and vector mesons, (2) charge radii of pion and kaon, (3) form factors of neutral and charged 
kaons, and (4) radiative decay widths for the V(P) — > P(V)j and P — > 77 transitions. These 
observables are calculated for the two potentials(HO and linear) to gauge the sensitivity of 
our results. 

Our calculation is carried out using the standard light-cone frame( q + = q° + q 3 = 0) with 
Qj_ = Q 2 = —q 2 - We think that this is a distinct advantage in the light-cone quark model 
because various form factor formulations are well established in the light-cone quantization 
method using this well-known Drell-Yan-West frame(g + = 0). The charge form factor of the 
pseudoscalar meson can be expressed for the '+ '-component of the current J^ 1 as follows 
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F(Q 2 ) = e q I(Q 2 , m g , m q ) + e q I(Q 2 , m q , m q ), 



(3.1) 



where e q (e q ) is the charge of quark(anti-quark) and 



'dk n u , , dk' ,, , „4 2 + k ± -k'_ 



I(Q 2 ,m q ,m q ) = dx dW^(s,kjW^*Qc,k', L ) - I ?— Mf> ( 3 - 2 ) 

7o 7 V «9X V «9X /^2 + k 2 7^2 + k /2 



with the definition of A and k'j_ given by 

A = xm q + (1 — x)m q , k'± = kj_ + (1 — x)q_i_. (3.3) 

Subsequently, the charge radius of the meson can be calculated by r 2 = — 6dF(Q 2 )/dQ 2 \Q2 =0 . 
Also, the decay constant fp of the pseudoscalar meson(P = ir, K) is given by 



fp = j/^j- 2 fdx [ d 2 kj^(x,k ± ) f (3.4) 
(2ti r 2 Jo J V dx yJA 2 + k 2 ± 

Since all other formulae for the physical observables such as the vector meson decay constants 
fv of V — > e + e~ , decay rates for the V(P) — > P(V)7 and P — > 77 transitions have already 
been given in our previous publication |J and also in other references(e.(?. Ref. 0), we do 
not list them here again. The readers are recommended to look at Refs. and M for the 
details of unlisted formulae. In Fig. 3, we show our numerical results of the pion form factor 
for HO (solid line) and linear (dotted line) cases and compare with the available experimental 
data [|3(]] up to Q 2 ~ 8 GeV 2 region. Since our model parameters of m u = 0.25 GeV and 
/3 U u= 0.3194 GeV for the HO case are same with the ones used in Refs. |7j] and [11 



our 



numerical result of the pion form factor is identical with the Fig. 2 (solid line) in Ref. [11 
In Figs. 4(a) and 4(b), we show our numerical results for the form factors of the charged 
and neutral kaons and compare with the results of vector model dominance (VMD) ||31|| . 
where a simple two-pole model of the kaon form factors was assumed, i.e., Fk+(k°)(Q 2 ) = 
e u(d) m2 j/( m t + Q 2 ) + e s m ^/( m2 i + Q 2 )- From Figs. 4(a) and 4(b), we can see that the neutral 
kaon form factors using the model parameters obtained from HO and linear potentials are 
not much different from each other even though the charged ones are somewhat different. 
The decay constants and charge radii of various pseudoscalar and vector mesons for the 



two potential cases are given in Table III and compared with experimental data [TT|,|3"2" 



While our optimal prediction of 5y was \5y\ = 4.2°[7.8°] for HOflinear] potential model, we 
displayed our results for the common 5y value with a small variation( i.e., \5y\ = 3.3° ± 1°) 
in Table III to show the sensitivity. The results for both potentials are not much different 
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from each other and both results are quite comparable with the experimental data. The 
decay widths of the V(P) — > P(V)7 transitions are also given for the two different potential 
models in Table IV. Although it is not easy to see which sign of 8y for HO potential model is 
more favorable to the experimental data, the positive sign of 5y looks a little better than the 
negative one for the processes of uj{4>) — > ?77 and r/ — > transitions. Especially, the overall 
predictions of HO potential model with the positive 5y seem to be in a good agreement 
with the experimental data. However, more observables should be compared with the data 
in order to give more definite answer for this sign issue of uj — <fi mixing angle. The overall 
predictions of linear potential model are also comparable with the experimental data even 
though the large variation of the mixing angle 5y should be taken into account in this case. 

In Table V, we show the results of P(= vr, 77,7/) — *> 77 decay widths obtained from 
our two potential models with the axial anomaly plus partial conservation of the axial 
current(PCAC) relations. The predictions of T](j]') — > 77 decay widths using PCAC are 
in a good agreement with the experimental data for both HO and linear potential models 
with 77 — 77' mixing angle, 6su(3) — —19°. The predictions of the decay constants for the 
octet and singlet mesons, i.e., r/ 8 and 770, are fs/ f n = 1.254[1.324] and fo/f n = 1.127[1.162] 
MeV for HO [linear] potential model, respectively. Our predictions of /§ and /o are not much 
different from the predictions of chiral perturbation theory [B3J reported as fs/f n — 1-25 and 
fo/fn — 1-04 ± 0.04, respectively. Another important mixing-independent quantity related 
to / 8 and f is the i?-ratio defined by 



R 



r(jj -> 77 ) r(rf -> 77) 



+ 



"'- KS. + S&). (3.5) 



r(vr^77) 3V/| /, 
Our predictions, R = 2.31 and 2.17 for HO and linear potential model cases, respectively, are 



quite comparable to the available experimental data |3^j3q| , -R ex p = 2.5±0.5(stat)±0.5(syst). 
Also, the Q 2 -dependent decay rates Tp^(Q 2 ) are calculated from the usual one-loop diagram 
^J^] and the results are shown in Figs. 5-7. Our results for both potential models are not only 
very similar with each other but also in a remarkably good agreement with the experimental 
data [36|-|3"8l up to Q 2 ~ 10 GeV 2 region. We think that the reason why our model is so 
successful for P — > 7*7 transition form factors is because the Q 2 -dependence(~ 1/Q 2 ) is due 
to the off-shell quark propagator in the one-loop diagram and there is no angular condition 
0] associated with the pseudoscalar meson. 
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IV. SUMMARY AND DISCUSSIONS 



In the light-cone quark model approach, we have investigated the mass spectra, mixing 
angles, and other physical observables of light pseudoscalar and vector mesons using QCD 
motivated potentials given by Eq.(2.2). The variational principle for the effective Hamilto- 
nian is crucial to find the optimum values of our model parameters. As shown in Figs. 1(a) 
and 1(b), we noticed that both central potentials in Eq.(2.2) are not only very similar to 
each other but also quite close to the ISGW2 ^ZB model potentials. In Figs. 1(a) and 1(b), 
we have also included the GI potential for comparison. Using the physical masses of (u, 0) 
and (rj, rj'), we were able to predict the u — <fi and 77 — 77' mixing angles as \8y\ ~ 4.2° [7.8°] 
and 6su(3) ~ —19.3° [—19.6°] for the HO [linear] potential model, respectively. We also have 
checked that the sensitivity of the mass spectra of (00, (p) to ~ 1°[5°] variation of 5y, i-e., 
from 5 V = 4.2°[7.8°] to 3.3° for HOflinear] potential case, is within 1%[5%] level. 

Then, we applied our models to compute the observables such as charge radii, decay 
constants, and radiative decays of P(V) — > V(P) r y* and P — > 77*. As summarized in Tables 
III, IV, and V, our numerical results for these observables in the two cases(HO and linear) 
are overall not much different from each other and are in a rather good agreement with 



the available experimental data [|Tj]]. Furthermore, our results of the i?-ratio presented in 
Eq.(3.5) are in a good agreement with the experimental data |34| , |3"5|| . The Q 2 dependence of 
P — > 77* processes were also compared with the experimental data up to Q 2 ~ 8 GeV 2 . The 
Q 2 -dependence for these processes is basically given by the off-shell quark propagator in the 
one-loop diagrams. As shown in Figs. 5-7, our results are in an excellent agreement with the 
experimental data pB|-P^]. Both the pion and kaon form factors were also predicted in Figs. 3 
and 4, respectively. We believe that the success of light-cone quark model hinges upon the 
suppression of complicated zero-mode contributions from the light-cone vacuum due to the 
rather large constituent quark masses. The well-established formulation of form factors in 
the Drell-Yan-West frame also plays an important role for our model to provide comparable 
result with the experimental data. Because of these successful applications of our variational 
effective Hamiltonian method, the extension to the heavy(6 and c quark sector) pseudoscalar 
and vector mesons and the ++ scalar mesons is currently under consideration. 

While there have been previous light-cone quark model results on the observables that 
we calculated in this work, they were based on the approach of modeling the wavefunction 
rather than modeling the potential. Our approach in this work attempting to fill the gap 



10 



between the model wavefunction and the QCD motivated potential has not yet been explored 
to cover as many observables as we did in this work. Nevertheless, it is not yet clear which 
sign of uj — mixing angle should be taken, even though the overall agreement between our 
HO potential model with the positive sign, i.e., 5y ~ 3.3° and the available experimental 
data seems to be quite good. If we were to choose the sign of X as X > in Eq.(B4), 
then the fact that the mass difference m u — m p is positive is correlated with the sign of 



the uj — 4> mixing angle [p9 |. In other words, > m p implies 8y > from Eqs.(B3)- 
(B5). Perhaps, the precision measurement of <fi —>■ 7/7 envisioned in the future at TJNAF 
experiment might be helpful to give more stringent test of Sy- In any case, more observables 
should be compared with the experimental data to give more definite assessment on this 
sign issue. 
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APPENDIX A: FIXATION OF THE MODEL PARAMETERS USING THE 

VARIATIONAL PRINCIPLE 



In this Appendix A, we discuss how to fix the parameters of our model, i.e., quark 
masses (m u , m s ), gaussian parameters( /3 uU = (5 U ^, j3 U g, f3 a g) and potential parameters (a, b, k) 
in V q q given by Eq.(2.2). In our potential model, the p — tt mass splitting is obtained by the 
hyperfine interaction, Vh yp . 

Our variational method first evaluates < ^[-Hq + Vq]\^/ > with a trial function <f) 1Q (k 2 ) 
that depends on the parameters (m, (3) and varies these parameters until the expectation 
value of H + Vq is a minumum. Once these model parameters are fixed, then, the mass 
eigenvalue of each meson is obtained by M qq =< V\[H + Vb]|* > + < ^\H hyp \^ > ft In 
this approach, we do not discriminate the gaussian parameter set (3 = (/3 m , j3 uS , j3 s s) by the 
spin structure of mesons. 

Let us now illustrate our detailed procedures of finding the optimized values of the model 
parameters using the variational principle: 

dp -°- (A) 

From Eqs.(2.1)-(2.2) and Eq.(Al), we obtain the following equations for HO and linear 
potentials: 

H.O. potent^ : b h = -[ ^=|, (A2) 

Linear potential : k = - | — 3^/?/' ^ ' 

where the subscript h[l] represents the HO [linear] potential parameters. Eqs.(A2) and (A3) 
imply that the variational principle reduces a degree of freedom in the parameter space. 
Thus, we have now four parameters, i.e., {m u , (3 U ^, a, ft(or b)}. However, in order to deter- 
mine these four parameters from the two experimental values of p and 7r masses, we need 
to choose two input parameters. These two parameters should be carefully chosen. Oth- 
erwise, even though the other two parameters are fixed by fitting the p and 7r masses, our 
predictions would be poor for other observables such as the ones in Sec. Ill as well as other 
mass spectra. From our trial and error type of analyses, we find that m u = 0.25 [0.22] GeV 



As we will see later, in our fitting of the p — tt splitting, the rather big mass shift due to the 
hyperfine interaction is attributed to the large QCD coupling constant, n = 0.3 ~ 0.6. 
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is the best input quark mass parameter for the HO [linear] potential among the widely used 
u(d) quark mass, m u = 0.22 GeV 0, 0.25 GeV §, and 0.33 GeV For the linear 

potential, the string tension b[ = 0.18 GeV 2 is well known from other quark model analyses 
P| j2^j25| commensurate with Regge phenomenology. Thus, we take m u = 0.22 GeV and 
hi = 0.18 GeV 2 as our input parameters for the linear potential case. However, for the HO 
potential, there is no well-known quantity corresponding to the string tension and thus we 
use the parameters of m u ^ = 0.25 GeV and (3 U ^ = 0.3194 GeV as our input parameters 
which turn out to be good values to describe various observables of both the it and p mesons 
for the gaussian radial wave function 

Using Eqs.(2.1),(A2) and (A3) with the input value sets of (l)(m u =0.25 GeV,/3 w -=0.3194 
GeV) for the HO potential and (2)(m u = 0.22 GeV, 6, = 0.18 GeV 2 ) for the linear potential, 
we obtain the following parameters from the p and it masses, viz., < ^/\H^j P ^\^/ >= M V ^ = 
m p(n)(P= Pseudoscalar and V= Vector): 



As shown in Fig. 1(a), it is interesting to note that our two central potentials, Coulomb plus 
HO (solid line) and Coulomb plus linear (dotted line) potentials, are not much different from 
each other and furthermore quite comparable to the Coulomb plus linear quark potential 
model suggested by Scora-Isgur(ISGW2) (long-dashed line for re = 0.3 and dot-dashed 
line for k = 0.6) up to the range of r < 2 fm. Those four potentials(H.O., Linear, and 
ISGW2) are also compared with the Godfrey and Isgur(GI) potential model |]J (short-dashed 
line) in Fig.l(a). The corresponding string tensions, i.e., /o( r ) = — dVoir) / dr , are also shown 
in Fig. 1(b). 

Next, among various sets of {m s1 (3 U s\ satisfying Eqs.(A2) and (A3), we find m s = 
0.48 [0.45] GeV and /? u5 =0.3419 [0.3886] GeV for HO [linear] potential by fitting optimally 
the masses of K* and K, i.e., M^s^ = ^\k*{k)- Once the set of {m s ,(3 uS } is fixed, then 
the parameters, (3 sS = 0.3681 [0.4128] GeV for HO [linear] potential, can be obtained from 
Eq.(A2)[(A3)]. Subsequently, and Af£ are predicted as 996 [952] MeV and 732 [734] 
MeV for HOflinear] potential, respectively. As shown in Fig.2(a)[2(b)], the solid, dotted and 
dot-dashed lines are fixed by the HO [linear] potential parameter sets of {m u , (3 u g}, {m s , j3 U g}, 
and /3 S s, respectively, and these three lines cross the same point in the space of b and k if 
the parameters in Table I are used. 



(1) H.O. potential: a h = -0.144 GeV, b h = 0.010 GeV 3 , K h = 0.607, 

(2) Linear potential : a t = -0.724 GeV, (3 ud = 0.3659 GeV, re, = 0.313. 



(A4) 
(A5) 
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We have also examined the sensitivity of our variational parameters and the correspond- 
ing mass spectra using a gaussian smearing function to weaken the singularity of 5 3 (r) in 
hyperfine interaction, viz., 



5 3 (r)^^exp(- -V 2 ). (A6) 

7T ' 

By adopting the well-known cut-off value of o = 1.8 Pl,|2"o| and repeating the same optimiza- 
tion procedure as the contact term(z.e., <5 3 (r)) case, we obtain the following parameters^ for 
each potential: 

H.O. potential : a h = -0.123 GeV, b h = 9.89 x 1(T 3 GeV 3 , n h = 0.636, (A7) 
Linear potential : a, = -0.7 GeV, k = 0.176 GeV 2 , k x = 0.332. (A8) 

The changes of other model parameters and mass spectra are given in Tables I and II. As 
one can see in Eqs.(A7)-(A8) and Tables I-II, the effects of smearing out S 3 (r) are quite 
small and the smearing effects are in fact negligible for our numerical analysis in Sec. III. 

APPENDIX B: MIXING ANGLES OF far/) AND 

In this appendix, we illustrate the mixing angles of (rj, rj') and (u, <p) by adopting the 
formulation to incorporate the quark-annihilation diagrams and the effect of SU(3) symmetry 
breaking in the meson mixing angles. 

The Eq.(2.5) satisfy the (mass) 2 eigenvalue equation 

M 2 \f t > =M 2 u \f l > (z = l,2). (Bl) 

Taking into account SU(3) symmetry breaking, we use the following parametrization for A4 2 
suggested by Scadron 



2 _(Ml- n + 2X V2XX \ 
M ~ { V2XX Ml + AX 2 ) ■ (B2) 

The parameter A characterizes the strength of the quark-annihilation graph which couples 
the 7=0 uu state to 7=0 uu, dd, ss states with equal strength in the exact SU(3) limit. 



For the sensitivity check of smearing out (5 3 (r)[Eq.(A6)], we kept /3 U(J = 0.3659 GeV for the linear 
potential case given by Eq.(A5) as an input value and checked how much bi changed. 
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The parameter X, however, pertains to SU(3) symmetry breaking such that the quark- 
annihilation graph factors into its flavor parts, with A , AX and XX 2 for the uu —>■ uu(dd), 
uu — > ss(or ss — > uu), and ss — > ss processes, respectively. Of course, X — > 1 in the SU(3) 
exact limit. Also, in Eq.(B2), M 2 n and M 2 S describe the masses of the corresponding mesons 
in the absence of mixing. 

Solving Eqs.(2.5),(Bl), and (B2), we obtain Eq.(2.6) and 

x = (Ml-Ml n ){Ml-Ml n ) 

x2 = 2(M% - M 2 S )(M 2 S - M 2 fi ) 

(M% - M^)(M 2 fl - M 2 n )> { ] 

2V2XX 

MS, + AJ-2A) - (B5) 
The Eqs.(B3) and (B4) are identical to the two constraints, Tr(.M 2 ) = Tr(M|.) and 
det(A4 2 )= det(M?.). The sign of 5 is fixed by the signs of the A and X from Eq.(B5). 
Also, since Eq.(B2) is decoupled from the subspace of (uu — dd)/\/2, the physical masses 
of m w and m p are confirmed to be the masses of M^ n and M^ n , respectively, as we used in 
Sec. II to fix the parameters (a, b, k). 

Given the fixed physical masses of M^ n = m n and M^ n = m p together with Mf. (i = 1,2), 
the magnitudes of mixing angles for 77 — 77' and u — (f) now depend only on the masses of M^ s 
and MY S , respectively, from Eq.(2.5). However, from Eqs.(B3)-(B5), one can see that the 
sign of mixing angle depends on the sign of parameter X. While Xp > is well supported 
by the particle data group |T7|](— 23° < ^sut^) ~ — 10°), the sign of X v is not yet definite 
at the present stage of phenomenology. Regarding on the sign of Xy, it is interesting to 
note that 5 V w -3.3°(= 9 SU{3) - 35.26°)(i.e., X v < 0) is favored in Refs. 0,^-^, while 
the conventional Gell-Mann-Okubo mass formula for the exact SU(3) limit (X — > 1) predicts 
5y 0° in the linear mass scheme and 5y ~ +3. 3° (i.e., Xy > 0) in the quadratic mass 
scheme ||17|| . Our predictions for the u — <p and 77 — 77' mixing angles are given in Sec. II. 

The corresponding results of the mixing parameters Xy(p) and Xym) in Eqs.(B3) and 
(B4) are obtained for the HO [linear] potential as follows 

X v = 0.57[0.73]m 2 GeV 2 , Xy = ±2.10[±3.08], 

X P = 13.5[13.3]m 2 GeV 2 , X P = 0.84[0.85]. (B6) 

Our values of Ay and Xp for both HO and linear potential cases are not much different from 
the predictions of Ref. pOfl . The reason why Xy is much smaller than Xp, i.e., Xp ~ 23[18]Ay 
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in our HO [linear] case and Xp ~ 18 Ay in Ref. ||20|| , may be attributed to the fact that in the 
quark- annihilation graph, the 1 annihilation graph involves one more gluon compared to 
the h annihilation graph. This also indicates the strong departure of rj — rj' from the ideal 
mixing. 
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TABLES 



TABLE I. Optimized quark masses {m q , m s ) and the gaussian parameters (5 for both harmonic 
oscillator and linear potentials obtained from the variational principle. q=u and d. 



Potential m q [GeV] 


m s [GeV] 


Ara[GeV] 


P sg [GeV] P qS [GeV] 


H.O. 0.25 


0.48 


0.3194 


0.3681 [0.3703] W 0.3419[0.3428] 


Linear 0.22 


0.45 


0.3659 


0.4128[0.4132] 0.3886[0.3887] 


[a] The values in parentheses are 


results from the 


smearing function Jl ,26 in Eq.(A6) 


instead of the contact term. 






TABLE II. Fit of the 


ground state 


meson masses with the parameters given in Table I. Un- 


derline masses are input data. The masses of {r),rf) and (co,<j)) were used to determine the mixing 


angles of rj — rf and u — <p 


respectively. 






1 Sq Experiment [Me V] 


H.O. 


Linear 3 Si 


Experiment H.O. Linear 


vr 135±0.00035 


135 


135 p 


770± 0.8 770 770 


K 494± 0.016 


470 [469] M 


478 [478] K* 


892± 0.24 875[875] 850[850] 


rj 547± 0.19 


547 


547 u 


782± 0.12 782 782 


i 958±0.14 


958 


958 <p 


1020±0.008 1020 1020 



^ The values in parentheses are results from the smearing function in Eq.(A6) instead 
of the contact term. 
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TABLE III. Decay constants and charge radii for various pseudoscalar and vector mesons. For 
comparison, we use \5y\ = 3.3° ± 1° for both potential cases. The experimental data are taken 
from Ref.[17], unless otherwise noted. 



5 V = -3.3° ±1° 5 V = +3.3° ± 1° 



Observables 


H.O. 


Linear 


H.O. 


Linear 


Experiment 


U [MeV] 


92.4 


91.8 


92.4 


91.8 


92.4±0.25 


fx [MeV] 


109.3 


114.1 


109.3 


114.1 


113.4± 1.1 


f P [MeV] 


151.9 


173.9 


151.9 


173.9 


152.8± 3.6 


fx* [MeV] 


157.6 


180.8 


157.6 


180.8 




f u [MeV] 


45.9±1.4 


52.6±1.6 


55.1±1.3 


63.1±1.5 


45.9± 0.7 


U [ MeV l 

r 2 [fm 2 ] 


82.6+ 0.8 


94.3+0.9 


76.7+ 1.0 


87.6+1.1 


79.1± 1.3 


0.449 


0.425 


0.449 


0.425 


0.432 ± 0.016 [32] 


r 2 K+ [fm 2 ] 


0.384 


0.354 


0.384 


0.354 


0.34± 0.05 [32] 


r 2 K0 [fm 2 ] 


-0.091 


-0.082 


-0.091 


-0.082 


-0.054± 0.101 [32] 



TABLE IV. Radiative decay widths for the V(P) — > P(V) 7 transitions. The mixing angles, 
@SU(3) = —19° for 77 — 77' and \Sy\ = 3.3° ± 1° for lj — <j>, are used for both potential models, 
respectively. The experimental data are taken from Ref.[17]. 



5 V = -3.3° ±1° S v = +3.3° ± 1 



Widths 


H.O. 


Linear 


H.O. 


Linear 


Experiment [keV] 


F(p± - 7r± 7 ) 


76 


69 


76 


69 


68 ±8 


r(a; — > 7r 7 ) 


730±1.3 


667±1.3 


730+1.3 


667=Fl-3 


717 ±51 


T(4> -> 7T7) 




5 1-2.6 
°- i +3.6 


5.6^§ 


K 1+3.6 

°- i -2.6 


5.8 ±0.6 


r(p -► 777) 


59 


54 


59 


54 


58 ± 10 


r(w -► 777) 


8.7+ 0.3 


7.9+0.3 


6.9+ 0.3 


6.3+0.3 


7.0 ± 1.8 


r(0 -► 777) 


38.7± 1.6 


37.8± 1.5 


49. 2± 1.6 


47.6± 1.5 


55.8 ±3.3 


r(V - P7) 


68 


62 


68 


62 


61 ±8 




4.9± 0.4 


4.5± 0.4 


7.6± 0.4 


7.0± 0.4 


6.1 ± 1.1 


r(0 -> 7/7) 


0.41+0.01 


0.39+0.01 


0.36+ 0.01 


0.34+0.01 


< 1.8 


r(K*° -► k° 7 ) 


124.5 


116.6 


124.5 


116.6 


117± 10 


r(is:*+ k+ 7 ) 


79.5 


71.4 


79.5 


71.4 


50 ± 5 



TABLE V. Radiative decay widths T(P — > 77) obtained by using the axial anomaly plus PC AC 
relations. 0su(3) = —19° for 77 — 77' mixing is used for both potential cases. The experimental data 
are taken from Ref.[17]. 



Widths 


H.O. 


Linear 


Experiment 


r(vr -► 77) 


7.73 


7.84 


7.8 ± 0.5[eV] 


-► 77) 


0.42 


0.42 


0.47 ± 0.05[keV] 


r(r?' - 77) 


4.1 


3.9 


4.3 ± 0.6[keV] 
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FIGURES 



Fig. 1(a). The central potential Vo(r) versus r. Our Coulomb plus harmonic oscilla- 
tor(solid line) and linear(dotted) potentials are compared with the quasi-relativistic poten- 
tial of ISGW2 [F2IJ model with k = 0.3(long-dashed line) and k = 0.6(dot-dashed line) and 
the relativized potential of GI ]]J model (short-dashed line). 

Fig. 1(b). The central force /o( r ) versus r. Our force for the linear potential is the 
same as that of ISGW2 with k = 0.3(dotted lines). The forces of GI and ISGW2 with 
k = 0.6 are the same with each other (dashed lines). Our force for the harmonic oscillator 
potential(solid line) are quite comparable with the other four forces up to the range of r < 2 
fm. 

Fig. 2 (a). The variational principle satisfying Eq.(A2). The solid, dotted, and dot- 
dashed lines are fixed by the sets of (m u ,P uU ), (m s ,/3 U s), and (m s ,(3 sS ), respectively. 

Fig.2(b). The variational principle satisfying Eq.(A3). The same line codes are used 
as Fig. 2 (a). 

Fig. 3. The charge form factor for the pion compared with data taken from Ref. 
PD| . The solid and dotted lines correspond to the results of harmonic oscillator and linear 
potential cases, respectively. 

Fig. 4 (a). Theoretical predictions of charged K + form factors using the parameters 
of both harmonic oscillator(solid) and linear(dotted) potentials compared with a simple 
two-pole VMD model (dot-dashed), F^f K0) = e u(d) ml / (ml + Q 2 ) + e.m^m^ + Q 2 ). 

Fig. 4(b). Theoretical predictions of neutral K° form factors. The same line codes 
are used as Fig.4(a). 

Fig. 5. The decay rate for the 7r — > 7*7 transition obtained from the one-loop diagram. 
Data are taken from Refs. [BE, 57]. 



Fig. 6. The decay rate for the 7] — > 7*7 transition obtained from the one-loop diagram. 
Data are taken from Refs. [|36|-|3"8|. 
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Fig. 7. The decay rate for the rf — > 7*7 transition obtained from the one- loop diagram. 
Data are taken from Refs. [^C -38 . 
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